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With strong gradients in the pedestal of high confinement mode (H-mode) fnsion plasmas, gy- 
rokinetic simulations are carried ont for the trapped electron and ion temperatnre gradient modes. 
A broad class of nnconventional mode structnres is fonnd to localize at arbitrary poloidal positions 
or with multiple peaks. It is found that these nnconventional ballooning structnres are associated 
with different eigen states for the most unstable mode. At weak gradient (low confinement mode 
or L-mode), the most unstable mode is usually in the ground eigen state, which corresponds to a 
conventional ballooning mode structure peaking in the outboard mid-plane of tokamaks. However, 
at strong gradient (H-mode), the most unstable mode is usually not the ground eigen state and the 
ballooning mode structure becomes unconventional. This result implies that the pedestal of H-mode 
could have better confinement than L-mode. 

PACS numbers: 52.35.Py, 52.30.Gz, 52.35.Kt 


Although numerous theoretical models have been 
suggested [T], a yet unexplained phenomenon in tokamak 
fusion plasmas, is the transition of low (L) to high (H) 
confinement states, where H-mode[2] has significant bet¬ 
ter confinement property than that of the L-mode. Un¬ 
derstanding of the H-mode physics is not only important 
to make controlled fusion more feasible, but also that 
the existence of and transitions among multi-equilibrium 
states are important fields of nonlinear physics in labo¬ 
ratory and the Universe. Drift wave turbulence is one 
of the major causes that leads to the anomalous trans¬ 
port widely observed in fusion and space plasmas O |1]. 
In order to control the turbulent transport, it is cru¬ 
cial to understand the underlying transport mechanism, 
which may vary for different types of instability that 
drive the turbulence. The correlation time and length 
are found to be closely related to the mode structure 
of the turbulence [S]. Therefore, the mode structure of 
the turbulence has a significant effect on the transport 
level [S]. 

In this Letter, we show that the linear properties of two 
major types of electrostatic micro-instabilities[3], namely 
the trapped electron mode (TEM) and ion temperature 
gradient (ITG) mode, are completely different in the 
H-mode (strong gradient) and L-mode (weak gradient) 
stages. With the conventional weak gradient, the mode 
structures for drift wave instabilities such as the ITG and 
TEM are of ballooning type, peaking at the outboard 
mid-plane of the tokamak (c.f., [Hi). This type of solu¬ 
tion has been intensively studied using the ballooning- 
representation[3 [TU] by reducing one two-dimensional 
(2D) real space eigen mode equation for the drift waves 
to two one-dimensional (ID) ballooning space eigen mode 
equations. For the 2D case we solve the eigen equation 
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FIG. 1: Typical electron temperature (also density) profiles 
nsed in this Letter. HL-2A L-mode edge plasma profile is 
weak gradient RoLf} < 40. HL-2A H-mode edge plasma 
profile is strong gradient RoLf} > 80. 


in the poloidal plane. For the ID case we solve the eigen 
equation in the parallel direction. The most unstable 
solutions in the ballooning space found in the past have 
usually the ballooning-angle parameter = om , which 
corresponds to the solution localized at the outside mid¬ 
plane, i.e., = 0 in our notation, where 0p is defined 

as the local peaking poloidal angle for the mode struc¬ 
ture. For this reason, many local eigenvalue codes such 
as HD7[T1] assume implicitly ifk = 0. The unconven¬ 
tional eigen modes with 9p ^ 0 have been recently discov¬ 
ered in the strong gradient parameter regime. Typically, 
|0p| ~ or < 7r/2 have been shown to exist [71151 lUl lU] . In 
this work, we find the most general unconventional eigen 
mode structures from first principle gyrokinetic simula¬ 
tions. The underlying physics is also explained and it has 
important implications for turbulent transport. 

We first obtain linear electrostatic results from global 
gyrokinetic particle simulation using the GTC code[T51 
nni with single toroidal mode number n. The simula¬ 
tion parameters and profiles are similar to that of the re¬ 
cent H-mode experiments of the HL-2A tokamak |171l24j : 
toroidal magnetic field Bq = l.SST, minor radius a = 
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FIG. 2: Conventional (a) and unconventional (b-i) 2D bal¬ 
looning structures of electrostatic potential in (X,Z) plane 
for TEM observed in GTC simulation, where (a) uses HL- 
2A tokamak edge weak gradient L-mode plasma parameter 
< 40) and (b)-(i) use edge strong gradient H-mode 
parameters {RqL^^ > 80). Collisions are only included in (e) 
and (g). 


40cm, major radius Rq = 165cm, safety factor q = 2.5 — 
3.0, magnetic shear s = 0.3— 1.0, Ro/Ln = 80 — 160 with 
Te{r) = Ti{r), and ne{r) = nj(r). = -{l/n){dn/dr) 
and Lt = —{l/T){dT/dr) are density and temperature 
gradient scale length. Typical electron temperature (also 
density) profiles used in this Letter are shown in Figj^ 
We start with rj = LnjLT = 1-0 for simplicity. Collisions 
are included in some cases but shown little influence to 
the general results. Under these parameters, no insta¬ 
bility or only weakly unstable mode can be found when 
the electrons are adiabatical. Thus, the major instability 
for these simulation parameters is the trapped electron 
mode. 

These TEM simulations show that both conventional 
and unconventional ballooning mode structures can ex¬ 
ist for various gradients and toroidal mode numbers 
(n = 5 — 30) , as shown by Fig[^ In these sub-figures, q 
profiles are similar. For Figl^-i, the global density (also 
temperature) profiles and toroidal mode numbers are not 
the same but all are under strong gradient. The novel 
features include: a), the mode can have anti-ballooning 
structure (i.e., \9p\ > 7r/2, e.g., Figj^); b). the mode 
can have multiple peaks (e.g., Figi^r Considering that 
the trapped particles are mainly located at the low mag¬ 
netic field side, i.e., the outboard side, the anti-ballooning 
structures of TEM are not expected. The 3D mode 





FIG. 3: The real part of Fourier 5<j)m.{r) for conventional and 
unconventional mode structures. The corresponding poloidal 
cross section mode structures of (a)-(d) (n = 20,10,5,10, 
respectively) are taken from Figi (a), (b), (g) and (i), re¬ 
spectively. The dash lines are corresponding rational surface 
positions r^, where nq{rs) = m. 


structure of the electrostatic potential can be represented 
by the Fourier series 6(j){r,9X) = e™*" 
where m is poloidal mode number. To explore the for¬ 
mation of these different eigenmode structures, we com¬ 
pute the S(pmi'r) for several typical conventional and un¬ 
conventional mode structures, as shown in Fig[^ For 
the conventional ballooning structure, the poloidal eigen 
modes 54>m{f) are almost radially symmetric (Gaussian- 
like) and positive in amplitude. And, dcfim has a large 
overlap with Sipm+i, i-e., 6(j)m — S4>m+i- However, for 
the unconventional structures, the poloidal eigen modes 
<5</'m('r) can be radially either symmetric or asymmetric, 
and the amplitude for each symmetric mode can be ei¬ 
ther positive or negative, as shown by Fig. 2b, c and d. 
Under stronger gradients the radial peaking position of 
dipniif) is also not at the corresponding rational surface 
position Tg any more, where nq(rs) = m. 

Next we consider ITG mode by reducing the den¬ 
sity gradient to Ro/Ln < 40 and keeping the other 
parameters the same as those for the TEM case, e.g., 
Rq/Lt > 80 and thus rji = L„/Lt > 2.0. To completely 
exclude the contribution of the kinetic electrons, we use 
adiabatic electrons in the simulations. It is found that 
the preceding unconventional mode structures still exist 
and exhibit even more structural variations. For exam¬ 
ple, the anti-ballooning structure is found for this ITG 
simulation, as is shown in Fig|^&b. Actually, the mode 
structure with global profiles and multi modes coexisting 
in the initial value simulation can be even more compli¬ 
cated. For example, two modes with similar growth rates 
can be excited in different radial locations, as shown in 
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FIG. 4: Unconventional ITG (n = 10) mode structures in 
GTG. (a & b) Anti-ballooning structure, (c & d) Two modes 
co-exist (or, one mode with two radius peaks) at different 
radius positions. One has 6p ~ 'k/2 and another has 6p ~ 
-tt/2. 


Fig|^&d. Multi modes coexist with close peaking po¬ 
sitions in the initial value simulation can also lead to 
ke., rotate poloidally with time. Thus, the 
unconventional mode structures are not limited to TEM 
and can be common for drift waves. 

These unconventional linear behaviors can be under¬ 
stood from the following eigenmode analysis. We start 
with the ITG eigen mode equation [51 [TU] 


2 ^ 

dx^ 




cos 9 - 


i sin 9 d 
ke dx 


cJ+r/s 


- 5(j){x,e) = 0 , 


( 1 ) 


where a = en/{qkePi), e„ = Lp/Ro, rjs = 1 + rji, 
X = r — Tg, Ts is the rational surface, w = Wr- -I- iy is 
the complex mode frequency normalized by the electron 
diamagnetic frequency, and the poloidal wave number 
ke = nq/r. Eq.([^ can be derived from the gyrokinetic 
theory with adiabatic electron assumption. The corre¬ 
sponding ID eigen mode equation in the ballooning space 
is 

{^£2 +k'^p'i[l + s'^{d-i9kf] + ^[cos9 
-hs(r?-■dfc)sini?]-h = 0, (2) 


where dk is the ballooning-angle parameter, which rep¬ 
resents an as yet undetermined radial wavenumber jlOj . 
The relation between the ballooning space electrostatic 
potential 6(i){'d,'dk) and real space d(j){x,9) can be found 
in Ref. [51 [rUj. Using the Fourier basis 6(j){x,6) = 


J2m ^4’me Eq.Q can be rewritten as the 2D eigen¬ 
mode equation 

klpfs^^-i^^ + £{z - mfHrn - ^ (1 - S-§£j5(l)m-l 
+ (^ +~ {zj^,+klp'£j54>m=(^, (3) 


where z = kgsx. To solve the eigenvalue problem of 
Eq.§, only a few number of m modes need to be kept 
for the solution to reach convergence. 

With suitable approximations (cf. Ref. [22]), both 
Eqs.([^ and (1^ can be reduced to the Weber equation 
u" + {J)x^ + a)u = 0 (here the argument x is d and z for 
Eqs.([^ and ^ respectively), which has solutions with 
the eigenvalues a(w) = i{2l + and eigenfunc¬ 

tions u{x) = Hi{iy/bx)e~''^^ where Hi is Uth Hermite 
polynomial and I = 0 , 1 , 2 ,..., which represent a series 
eigenstates. With the original equations, i.e., Eqs. (§ and 
(§, which can only be solved numerically, the eigenstates 
take a more complicated form. 

Eqs. (§ and (|^ can be solved numerically by trans¬ 
forming it to a matrix eigenvalue problem as uj^M^X -|- 
w^M 2 X+u}MiX+MqX = 0, with X is the discrete rep¬ 
resentation of the electrostatic potential. We use finite 
difference to discretize the system, which yields sparse 
matrices for Mi {i = 0,1,2,3). Using the companion 
matrix method, the nonlinear eigenvalue problem can be 
transformed to a standard eigenvalue problem as AY = 
wBY, where Y = [Xi,X 2 ,X^\ = [X,ujX,w‘^XI 
A = [0,1,0-0,0,1,-Mo,-Mr,-M 2 I B = 

[/, 0,0; 0,1,0; 0,0, M 3 ], and I and O are unit and 
null matrix respectively. Thus all the solutions of this 
eigen system can be obtained (c.f., [T5] for details of sim¬ 
ilar treatment). The advantage of this method is that 
it can show the complete set solutions of the discrete 
eigen system and help us to understand the distribu¬ 
tion of eigenvalues in the complex plane. The solution 
in Refs. 130 using iterative solver is actually just one of 
the solutions obtained here and may not be the most un¬ 
stable or most important, which depends heavily on the 
initial guess. This companion matrix method has been 
verified by comparing the numerical solutions with that 
from the shooting method and the analytical solution for 
the Weber equation. 

By solving Eq. (U in the ID ballooning space, the un¬ 
conventional ballooning structures occur when either the 
most unstable solution is not the ground eigen state 
(Z ^ 0), or the ballooning angle dk ^ 0. Both of these 
two conditions can be met in the strong gradient regime. 
The most unstable solution with 7 ^ 0 has been dis¬ 
cussed by others (c.f., [13 US]). Here we focus on the 
unconventional ballooning structure caused by the non¬ 
ground eigen state. The following parameters are used to 
solve Eq.([^: s = 0.8, kgpi = 0.4, q = 1.0, rjg = 3.0 and 
dfc = 0. As is known from the aforementioned analytical 
analogy, Fig[5] shows that a series of solutions can exist 
for Eq. where R/Lj' and R/L^ are changed simulta¬ 
neously to ensure rj = L^ILt = const., and one should 
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FIG. 5: In Eq. series of solutions exist. For weak gradient 
(£„ = 0.5), the most unstable solution (red ’x’) is the ground 
state (a&b), which is the conventional ballooning structure. 
For strong gradient (e„ = 0.2), the most unstable solution 
(red ’x’) is not the ground state (c&d), which represents the 
unconventional ballooning structure. 




FIG. 6: Typical unconventional mode structures from 2D 
eigen solution for Eq.([3l. (b) is similar to Fig[^c&d), and 
(c&d) are similar to Fig ^d&c) 


also note that the frequency is normalized by the elec¬ 
tron diamagnetic frequency. For the weak gradient case 
(e„ = 0.5), we find that the most unstable solution is 
the ground state (Figl^), which is the conventional bal¬ 
looning structure (Figj^). For the strong gradient case 
(e„ = 0.2), the most unstable solution is not the ground 
state (Figj^&d), which corresponds to the unconven¬ 
tional ballooning structure. More detailed analysis[2T] of 
Eq.ll]) for present discussion of the unconventional mode 
structure can be obtained by extension of Refs. uni^. 


We have demonstrated that, with strong gradient 
the most unstable solution can shift from ground state 
to other non-ground states, which is analogous to the 
quantum jump between energy levels. Physically, the 
jump behavior can be understood from the effective 
potential [20]. The jump happens from one potential well 
to another, which leads to different energy levels. It is 
not transparent that the non-ground eigen state in the 
ID ballooning space corresponds to the unconventional 
mode structure in the 2D poloidal plane. Next we con¬ 
firm this link by showing that the non-ground 2D eigen 
state solved from Eq.([^ can form the unconventional 
mode structures observed in the preceding gyrokinetic 
simulation. The solutions in Refs. [aiiiiii] are just a weak 
asymmetric solutions of our series solutions. Almost all 
the mode structures in Figsand have also been found 
in the 2D eigen solutions of Eq.([^. Two examples are 
shown in Figj^ Therefore, conventional and unconven¬ 
tional series solutions have been found in both 2D eigen 
solver and GTC initial simulations. The condition for 
the jump of the most unstable eigen state to non-ground 
state is e„ < Cc, where Cc is a critical gradient parameter 
which depends on other parameters. In GTC simulations 
of the HL-2A parameters, the typical critical density (or 
temperature) gradient value is Ro/Ln = 40 — 120. 

The results from the gyrokinetic simulation and eigen 
mode analysis show that the unconventional mode struc¬ 
tures exist mainly in the strong gradient regime or the 
H-mode. In the weak gradient regime or L-mode, con¬ 
ventional mode structures still prevail. This can indi¬ 
cate different transport behavior between H-mode and 
L-mode[23]. In the conventional ballooning structure, 
the neighboring Fourier modes 64>m — 54>m+i, the effec¬ 
tive correlation length may be estimated as the width of 
radial envelope of the modes, say, AA. Whereas, in the 
unconventional ballooning structures, especially for anti¬ 
ballooning structure, 54>m — —54>m+i can occur, i.e., a 
180° phase shift for the neighboring Fourier modes, which 
can change the effective correlation length to the distance 
of neighboring mode-rational surfaces Arg. Gonsidering 
that Arg <C AA, we can expect that the H-mode can 
have better confinement. 

To summarize, a broad class of unconventional balloon¬ 
ing modes are found for electrostatic drift waves (TEM 
and ITG) by the gyrokinetic simulation, which is shown 
to be common in the strong gradient regime. These un¬ 
conventional mode structures are shown to correspond to 
the non-ground-state solutions of the eigen mode equa¬ 
tion. These results may have important implications 
for the turbulent transport in tokamaks, i.e., the turbu¬ 
lent transport mechanism in the H-mode can be rather 
different from that in the L-mode, which requires fur¬ 
ther investigation by self-consistent nonlinear gyrokinetic 
simulations [23] . 
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